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Abstract. One of the main challenges of modern cosmology is to investigate the nature of 
dark energy in our Universe. The properties of such a component are normally summarised 
as a perfect fluid with a (potentially) time-dependent equation-of-state parameter 'w{z). We 
investigate the evolution of this parameter with redshift by performing a Bayesian analysis 
of current cosmological observations. We model the temporal evolution as piecewise linear 
in redshift between 'nodes', whose tf- values and redshifts are allowed to vary. The optimal 
number of nodes is chosen by the Bayesian evidence. In this way, we can both determine 
the complexity supported by current data and locate any features present in w{z). We 
compare this node-based reconstruction with some previously well-studied parameterisations: 
the Chevallier-Polarski-Linder (CPL), the Jassal-Bagla-Padmanabhan (JBP) and the Felice- 
Nesseris-Tsujikawa (FNT). By comparing the Bayesian evidence for all of these models we 
find an indication towards possible time-dependence in the dark energy equation-of-state. It 
is also worth noting that the CPL and JBP models are strongly disfavoured, whilst the FNT 
is just significantly disfavoured, when compared to a simple cosmological constant w = —1. 
We find that our node-based reconstruction model is slightly disfavoured with respect to the 
ACDM model. 

Keywords: Equation of State, Dark Energy, Cosmological Parameters from CMBR, Bayesian 
Analysis 
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1 Introduction 

Over the past decade, one of the most pressing goals of modern cosmology has been to explain 
the accelerated expansion of the Universe [1, 2]. Considerable observational and theoretical 
effort has been focused on understanding this remarkable phenomenon. It is often postulated 
that an exotic new source of stress-energy with negative pressure may be responsible for the 
cosmic acceleration: such a component is called dark energy (DE). 

The dynamical properties of dark energy are normally summarised as a perfect fluid with 
(in general) a time-dependent equation-of-state parameter w{z), defined as the ratio of its 
pressure to its energy density. The simplest proposal, namely a cosmological constant A, is 
described by the redshift independent w = —1. Alternative cosmological models that de- 
viate from standard ACDM, but still lead to an accelerating Universe, include: K-essence, 
quintessence and non-minimally coupled scalar fields [3-6], braneworld models [7], modified 
gravity [8-12], interacting dark energy [13-15], anisotropic universes [16-18], amongst many 
others [19-26]. In the absence of a fundamental and well-defined theory of dark energy, 
w{z) has been parameterised in a number of different ways, including: the CPL, JBP and 
FNT models [27-30], the Hannestad and Wetterich parameterisations [31, 32], polynomial, 
logarithmic and oscillatory expansions [33-35], Kink models [36], and quite a few others 
[37]. The a priori assumption of a specific model or the use of particular parameterisa- 
tions can, however, lead to misleading results regarding the properties of the dark energy. 
Hence, some studies instead perform a direct, model-independent ('free-form') reconstruction 
of w{z) from observational data, using, for instance, a principal component analysis [38-43], 
maximum entropy techniques [44], binning w{z) in redshift space [45, 46], non-parametric 
approaches [47-52] and several other techniques [53-67] 

In this paper we explore the possible dynamical behaviour of the dark energy based 
on the most minimal a priori assumptions. Given current cosmological observations and 
using the Bayesian evidence as an implementation of Occam's razor, we select the preferred 
shape of w{z). Our method considers possible deviations from the cosmological constant 
by modelling w{z) as a linear interpolation between a set of 'nodes' with varying w-values 
and redshifts (in the most general case). An advantage of this method is that the number 
of nodes is directly chosen by the model Bayesian evidence. This reconstruction process is 
essentially identical to the approach used previously to recover the preferred shape of the 
primordial spectrum of curvature perturbations P{k) [68]. For comparison, we also consider 
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some existing models that propose a parameterised functional form for w{z), namely the 
CPL, JBP and FNT models. For each model we compute its evidence and, according to the 
Jeffreys guidelines, we select the best model preferred by the data. 

The paper is organised as follows: in the next Section we describe the data sets and 
cosmological parameters used in the analysis. We then describe the form of existing param- 
eterisations used by other authors and define the reconstruction used in this paper. The 
resulting parameter constraints and evidences for each model are then discussed. Finally, in 
Section 3, based on Jeffrey's guidelines, we decide which model provides the best description 
for current observational data and present our conclusions. 

2 Analysis 

The data-sets considered throughout our analysis include temperature and polarisation mea- 
surements from the 7- year data release of the Wilkinson Microwave Anisotropy Probe (WMAP; 
[69]), together with the 148 GHz measurements from the Atacama Cosmology Telescope 
(ACT; [70]). In addition to CMB data, we include distance measurements of 557 Super- 
novae la from the Supernova Cosmology Project Union 2 compilation (SCP; [71]). We also 
incorporate Baryon Acoustic Oscillation (BAO; [72]) measurements of distance, and baryon 
density information from Big Bang Nucleosyntesis (BBN; [73]), and impose a Gaussian prior 
using measurements of the Hubble parameter today Hq, from the Hubble Space Telescope 
key project (HST; [74]). 

We consider purely Gaussian adiabatic scalar perturbations and neglect tensor contri- 
butions. We assume a flat CDM universe^ described by the following parameters: Ob/i^ and 
^DM^^ are the physical baryon and dark matter densities, respectively, relative to the critical 
density (h is the dimensionless Hubble parameter such that Hq = lOO/i kms^^Mpc^^), 6 is 
100 X the ratio of the sound horizon to angular diameter distance at last scattering surface, r 
is the optical depth at reionisation, Ag and Ug are the amplitude of the primordial spectrum 
and the spectral index respectively, measured at the pivot scale ko = 0.002 Mpc~^. Aside 
from the Sunyaev-Zel'dovich (SZ) amplitude Asz used by WMAP analyses, the 148 GHz 
ACT likelihood incorporates two additional secondary parameters: the total Poisson power 
Ap at I = 3000 and the amplitude of the clustered power Ac. To describe the overall shape of 
the dark energy equation-of-state parameter w{z) in our nodal reconstruction, we introduce 
a set of amplitudes Wz^ at determined positions Zj. The CPL and JBP models each depend 
upon just two parameters: wq and Wa, whereas the FNT parameterisation depends upon 
four parameters: wq, Wa, t and at- The assumed flat priors on the parameters of each w{z) 
reconstruction are discussed below. 

To carry out the exploration of the parameter space, we input w{z) into a modified 
version of the C AMB code [79] , which implements a parameterised post-Friedmann (PPF) 
presciption for the dark energy perturbations [80]. Then, we incorporate into the COSMOMC 
package [81] a substantially improved and fully-parallelized version of the nested sampling 
algorithm MultiNest [82, 83]. The MultiNest algorithm increases the samphng efficiency 
for calculating the evidence and allows one to obtain posterior samples even from distributions 
with multiple modes and/or pronounced degeneracies between parameters. The Bayes factor 

^The possibility of a dynamical dark energy in a curved universe has also been considered by, i.e. [75-78]. 
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Bij , or equivalently the difference in log evidences In Zi — In Zj , provides a measure of how 
well model i fits the data compared to model j [84-87]. A suitable guideline for making 
qualitative conclusions has been laid out by Jeffreys [88]: if Bij < 1 model i should not be 
favoured over model j, 1 < Bij < 2.5 constitutes significant evidence, 2.5 < Bij < 5 is strong 
evidence, while Bij > 5 would be considered decisive. 

2.1 Nodal Reconstruction I 

We first perform the reconstruction of 'w{z) by parameterising it as piecewise linear between 
a set of nodes with variable amplitudes (t(72.-values), but with fixed, equally-spaced redshifts. 
Throughout, we bear in mind that current relevant information, mainly coming from SN la, 
is encompassed between the present epoch Zmm = and Zmax = 2. At higher redshifts there 
is no substantial information to place strong constraints on dark energy, thus beyond z > 2 
we assume "^(2;) to be constant, with a value equal to that at Zmax- At each node, we allow 
variations in amplitudes w^. with a conservative prior w^. = [—2, 0]. Our description of w{z) 
can be summarised as: 

z = 

z e {zi} (2.1) 
z>2 

and with linear interpolation for < Zi < z < Zj+i < 2. 

While the use of linear interpolation between nodes may seem crude, we have shown in a 
previous work [68] that the use of smoothed interpolation functions, such as cubic splines, 
can lead to significant spurious features in the reconstruction, thus leading to poor fits to 
observational data and also unrepresentative errors. 

We perform all of our model comparisons with respect to the simplest explanation of 
dark energy, namely a cosmological constant, which is specified by a redshift-independent 
w = —1. First, we consider deviations of the ACDM model by letting the equation-of- 
state parameter vary only in amplitude: 'w{z) = wq = constant (see Figure 1(a)). The 
incorporation of two or more parameters, as in models (b) and (c) respectively, allows us to 
test the dark energy time-evolution. Figure 1 also includes the ID marginalised posterior 
distribution for the corresponding amplitude at each node and for each reconstruction. In the 
top label of each model we have included the Bayes factor compared to the ACDM model. 

In model (b), we notice the overall shape of w{z) includes a slight positive tilt and a 
narrow waist located at z ~ 0.3. It is also observed that at the present epoch w{z = 0) < —1 
is slightly favoured, while at higher redshifts 'w{z) > —1 is preferred, hence, the reconstructed 
w{z) exhibit a crossing of the line w = —1. The crossing of the phantom divide line w = —1 
(PDL) , plays a key role in identifying the correct dark energy model [89] . If future surveys 
confirm its existence, single scalar field theories (with minimal assumptions) might be in 
serious problems as they cannot reproduce this essential feature, and therefore alternative 
models should be considered, e.g. scalar-tensor theories [90, 91], braneworld models [92, 93], 
f{R) gravity [10-12, 94]. To continue with our reconstruction process, we then place a third 
point (c) midway between the two existing nodes in (b). This model mimics a running 
behaviour by allowing slight variations in the interpolated slopes between the three nodes. 
The freedom in its shape, together with the very weak constraints at high redshifts, lead to 
a w{z) with slight negative slope at early times, in contrast to model (b). Furthermore the 
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(a) Bi^A = -2.19 ±0.35 




z 



(b) B2,A = -2.34 ± 0.35 




(d) Bi^A = -1.57 ±0.35 




Figure 1. Left: Reconstruction of the dark energy equation-of-state parameter modelled as piecewise 
linear between nodes that may vary in amplitude Wi but are fixed in redshift z, showing the mean 
amplitude values and their corresponding la error bands. The colour-code shows In(likelihood), 
where lighter regions represents an improved fit. Right: ID marginalised posterior distribution of the 
amplitudes Wi at each z-node (shown in the right-top corner), in each reconstruction. The top label 
in each panel denotes the associated Bayes factor respect to the ACDM model. 
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presence of a small bump in the resulting w{z) at z ~ 1 (see Figure 1 (c)) might point to 
some weak departure from the cosmological constant w = —1. 

We can continue this process of adding more nodes but always using the Bayesian evidence 
to penalise any unnecessary inclusion of model parameters. The inclusion of a fourth stage 
with z-space split into three equally spaced regions is given by model (d). At low redshifts 
the shape of the equation of state is well constrained with tight error bands on each node, 
whereas at high redshifts the error bands again indicate the lack of sufficient data to provide 
strong constraints. Notice also the increased error bands due to the addition of further 
nodes and (anti-)correlations created between them: for instance, the posterior distribution 
of the amplitude Wz^ at z = is broadened as the number of nodes is increased. At this 
stage, the evidence has flattened off, and so it seems reasonable to stop adding parameters 
in the reconstruction process at this point. The constraints on the w^. amplitudes used on 
each reconstruction are given by (for two-tailed distributions 68% C.L. are shown, whilst for 
one-tailed distributions the upper 95% C.L.): 

(a) ?i;o = -1.02 ±0.07, 

(6) tt;^=o = -1-09 ±0.14, tf^>2 = -0.83 ± 0.39, 

(c) Wz=o = -l.U±0.17, u'^=i = -0.73 ± 0.33, Wz>2 < -0-Q5, 

(d) Wz=o = -1.18 ± 0.20, Wz=om = -0.78 ± 0.30, ^^=1.33 = 1-03 ± 0.53, Wz>2 < -0.62. 

The models used in the reconstruction of w{z) are assessed according to the Jeffreys 
guideline. The Bayes factor between the ACDM model and the one- node model Bi^a = 
—2.19 ±0.35 points out that ■w{z) = wq (where wq is a free constant), is strongly disfavoured 
when compared to the cosmological constant, similarly, when two independent nodes are 
used B2,A = —2.34 ± 0.35. Thus, parameterisations that contain one or two parameters 
are not able to provide an adequate description of the behaviour of w{z), and hence are 
strongly disfavoured by current observations. The addition of nodes in the third and fourth 
stage provides more flexibility in the shape of the reconstructed w{z). Thus, the evidence for 
these models shows an improvement, compared to the first and second models, indicating the 
possible presence of some features in the time evolution of the equation-of-state parameter. 
Nonetheless, when they are compared to ACDM, they are still marginally disfavoured: B^^a = 
-1.70 ± 0.35 and S4,A = -1-57 ± 0.35. 

2.2 Nodal reconstruction II 

We previously reconstructed w{z) by placing nodes at particular fixed positions in z-space. 
However, to localise features, we now extend the analysis by also allowing the z-position of 
each node to move freely. In particular, we again fix two z-nodes at sufficiently separated 
positions Zmin = and Zmax = 2, but now place inside additional 'nodes' with the freedom 
to move around in both position Zi and amplitude Wz^ ■ This method has the advantage that 
we do not have to specify the number and location of nodes describing w{z); indeed, the 
form of any deviation from flat w{z) can be mimicked through a change in the amplitudes 
and/or positions of the internal nodes. Also, the reduced number of internal nodes avoids 
the creation of wiggles due to high (anti-)correlation between nodes, which might lead to a 
misleading shape for w{z). We use the same priors for the amplitudes Wz^ = [—2,0] as we 
adopted in Section 2.1. Hence, for this type of nodal-reconstruction the equation of state is 
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-1.27 ±0.35 




Figure 2. Left: Reconstruction of the dark energy equation-of-state parameter w{z) using one- 
internal-node (top) and two-internal z-nodes (bottom) that move freely in both amplitude Wi and 
redshift z^. Right: corresponds to the ID and 2D marginalised posterior distribution of the amplitudes 
and z-node positions in each reconstruction. The colour-code indicates the In(Likelihood), where 
lighter regions represents an improved fit, and the top label in each panel denotes the associated 
Bayes factor with respect to the ACDM model. 



described by 



w 



z = 

0<Zi< Zi+i < 2 
z> 2 



(2.2) 



and with linear interpolation for < zi < Zj+i < 2. 

Figure 2 illustrates the reconstruction of w{z) from the mean posterior estimates for each 
node, with la error bands on the amplitudes (left). Also plotted are the ID and 2D 
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marginalised posterior distributions on the parameters used to describe w{z) (right). The 
reconstructed shape for the two-internal-node model (middle panel) resembles the form ob- 
tained in Figure 1(c), but now with a turn-over shifted to earlier times. A similar turn-over 
has been found using principal component analysis by [40, 41]. The narrow waist at z ~ 0.3 
is also noticeable, where the SNe constraints seem to be tightest. For the one and three- 
internal-nodes case (top and bottom panel of Figure 2), we observe w{z) has essentially the 
same behaviour as in the two-internal-node model, being the preferred model. Finally, a 
common feature throughout all the reconstructed equation of state w{z) is observed: the 
presence of the crossing PDL within the range < z < 0.5. The constraints on the Wz^ 
amplitudes used on each reconstruction are given by (for two-tailed distributions 68% C.L. 
are shown, whilst for one-tailed distributions the upper 95% C.L.): 

(zi) w,=o = -1.14 ± 0.18, wo<z<2 > -1-39 ± 0.35, w,>2 < -0.70, 

{z2) Wz=o = -1.18 ± 0.26, tt;o<^<i = -0.83 ± 0.29, wi<z<2 = 1-02 ± 0.52, Wz>2 < -0.63, 
(zs) w,=o = -1-07 ± 0.36, ?i;o<.<o.66 = -0.98 ± 0.29, u;o.66<^<i.33 = -0.84 ± 0.47, 
w'i.33<.<2 = -1.02 ± 0.55, Wz>2 < 0.63. 

The similar shape of the three models are in good agreement with their Bayes factor: 
i3z2,zi = +0.46 ± 0.35, 13z3,z2 = —0.14 it 0.35. According to the Jeffreys guideline, even 
though the two internal-node model contains more parameters, it is significantly preferred 
over the models with one and two fixed-nodes, i.e. Bz2,2 = +1.53 ± 0.35. However, when 
compared to the cosmological constant model the Bayes factor is too small to draw any 
decisive conclusions: Bz2,A = —0.81 ± 0.35. Thus we conclude that the internal-node models 
might be considered as viable models to characterise the dark energy dynamics. As seen 
in Figure 2, the Bayesian evidence has reached a plateau and thus we cease the addition of 
further nodes. 

2.3 CPL and JBP parameterisations 

In this section we examine some existing parameterised models for w{z) and compare these 
to our nodal reconstructions. In particular, we consider the simple parameterised description 
introduced by Chevallier-Polarski-Linder (CPL; [27, 28]), that has the functional form: 

w{z) =wo + Wa——, (2.3) 
1 -\- z 

where the parameters wq and Wa are real numbers such that at the present epoch w\z=o = wq 
and dw/dz\z=o = —Wa, and as we go back in time w{z » 1) ~ tt;o + Wa- Thus, we limit the 
CPL parameters by the flat priors wq = [—2,0] and Wa = [—3,2]. 

We also consider the parameterisation suggested by Jassal-Bagla-Padmanabhan (JBP; [29]): 

w{z) = Wo + Waj—^ — ^- (2.4) 

(1 + zy 

In this model, the parameter wq determines the properties of w{z) at both low and high 
redshifts: 'w{z = 0) = wo and w{z S> 1) ~ wq. To explore the parameter space we consider 
the following flat priors on the JBP parameters: wq = [—2, 0] and Wa = [—6, 6]. 
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(CPL) ScPL,A = -2.84 ± 0.35 




(JBP) BjBP,A = -2.82 ± 0.35 




Figure 3. Reconstruction of the dark energy equation of state 'w{z) assuming tire Chevallier-Polarski- 
Linder (top) and the Jassal-Bagla-Padmanabhan parameterisation (bottom), along with their corre- 
sponding 2D constraints with la and 2a confidence contours (right panel). The colour-code indicates 
the In(Likelihood), where lighter regions represents an improved fit; the top label in each panel de- 
notes the associated Bayes factor with respect to the ACDM model. Dotted lines indicate the priors 
choice. 

Figure 3 shows 2D joint constraints, with la and 2a confidence contours, for the parame- 
ters used to describe the CPL and JBP models, and the resulting shape of w{z) corresponding 
to the mean posterior estimates of wo and Wa- In each panel we have included the Bayes 
factor compared to the ACDM model. Both of the models are in good agreement with a 
simple cosmological constant. The current constraints for the CPL and JBP parameters are 
essentially as we expected: 

(CPL) t(;o = -1.11 ±0.17, = 0.34 ± 0.60, 
(JBP) Wo = -1.21 ± 0.26, Wa = 1.28 ± 1.62. 

Given that the CPL and JBP parametererisations depend upon just two parameters, 
they seem to not posses enough freedom to capture local features of w{z), i.e. the CPL 
model does not exhibit a turn-over, see Figure 3. This is reflected in the large difference 
in the Bayesian evidence for this model compared to that of the cosmological constant: 
■ScPL.A = -2.84 ± 0.35 and ^jbp,a = -2.82 ± 0.35. Li fact, the CPL equation of state 
looks similar to that obtained in Figure 1 (b), confirming our results. An important point to 
emphasise is that, for the chosen priors, Bcpl,z2 = —2.03 ± 0.35 and i3jBP,^2 = —2.01 ± 0.35, 
indicating that both models are strongly disfavoured in comparison to the internal-node re- 
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Table 1. Robustness of the CPL and JBP models over small variations of the prior range. The 
associated Bayes factor in each model is compared with respect to the ACDM model. 



Prior 


^CPL.A 


Prior 


^JBP,A 


Wo,Wa 




Wo,Wa 




[-1.5,-0.7], [-3,2] 


-1.84 ±0.35 


[-1.8,-0.6], [-6,6] 


-2.35 ±0.35 


[2,0], [-0.5,1] 


-2.11 ±0.35 


[-2,0], [-1,4] 


-1.82 ±0.35 


[-1.5,-0.7], [-0.5,1] 


-1.39 ±0.35 


[-1.8,-0.6], [-1,4] 


-1.51 ±0.35 


[-1.3,-1], [0,1] 


-0.26 ±0.35 


[-1.4,-1.1], [0,3] 


-0.54 ±0.35 



construction, shown in Figure 2. 

To illustrate the robustness of the model to small variations of the prior range, we com- 
pute the Bayesian evidence using different sets of priors, shown in Table 1; the prior ranges 
are illustrated with dotted lines in Figure 3. The reader will observe that even though the 
priors, in the first three choices, have been shrunk to within the region of the 2a contours, 
the Bayes factor still disfavours significantly both the CPL and JBP parameterisations com- 
pared to the cosmological constant and the two-internal-node reconstruction. With respect 
to the extremely small prior (last row of Table 1), we notice that the JBP model does not 
contain the cosmological constant wq = —1. Its Bayes factor compared to the ACDM model 
•SjBP.A = —0.54 ± 0.35, shows that models with w{z = 0) < —1.1 might provide a good 
description for the current state of the Universe. 

2.4 FNT parameterisation 

We have observed that two-parameter functions are not, in general, sufficient to recover the 
evolution of the dark energy w{z), obtained previously in the reconstruction process. As an 
alternative to the CPL and JBP functional form, we consider a more general parameterisation 
introduced by Felice-Nesseris-Tsujikawa (FNT, [30]), which allows fast transitions for the dark 
energy equation of state: 

w{a) = Wa^{wQ-Wa) , ^ ^ , (2.5) 

where a = 1/(1 + z), > and r > 0. The parameter wq determines the w{a) properties at 
present time wi^ = w{a = 1), whereas Wa the asymptotic past Wa = w{a <C 1). In this model, 
the equation of state w{a) has an extremum at a^, = at with value 

w{a^)=Wp + - _^ . (2.6) 

Based on the assumptions given by [30], we explore the cosmological parameter-space using 
the following flat priors: wq = [—2,0], Wa = [—2,0], at = [0, 1] and r = [0, 1], using a full 
Monte-Carlo exploration. We leave the analysis of the robustness of this model under small 
variations on the priors, for a future work. 
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(FNT) ^FNT,A = 



-1.68 ±0.35 




















1^ 






1/.)/ 



Figure 4. Reconstruction of the dark energy equation of state w{z) assuming the Fehce-Nesseris- 
Tsujikawa paramterisation (left panel), along with their corresponding ID, and 2D constraints with 
la and 2a confidence contours (right panel). The colour-code indicates the In(Likelihood), where 
lighter regions represents an improved fit; the top label in the panel denotes the associated Bayes 
factor with respect to the ACDM model. 



In Figure 4 we plot 2D joint constraints, with la and 2a confidence contours, for the 
parameters used to describe the FNT model, and its corresponding reconstruction of w{z). 
We observe that the FNT model is in good agreement with a simple cosmological constant 
w{z) = —1, with current constraints: 



(FNT) 1/70 = -1.19 ±0.32, = -0.94 ± 0.15. 

Given that the best-fit values of wq and Wa are very similar, the second term on the left hand 
side of (2.6) is almost negligible. This results in essentially unconstrained values for at and r, 
and so Wa becomes the dominant term in the dynamics of w{z). We have found that the FNT 
model shares a similar feature common throughout all the models: w{z = 0) < w{z ^ 1), in 
agreement with our previous results. The best-fit form of w{z) presents a maximum value 
given by u)(a*) = —0.95 located at z* = 1/a* — 1 = 1.59. On the other hand, the top label 
of Figure 4 shows the Bayes factor compared to the ACDM model: i3FNT,A = —1.68 ± 0.35. 
That is, the FNT model improves on the Evidence computed from the CPL and JBP models, 
however the inclusion of twice the number of parameters makes it significantly disfavored 
when compared to the cosmological constant w{z) = —1, and indistinguisable compared to 
our node-base reconstruction, i.e. i3FNT,z2 — —0-82 ± 0.35. 

3 Discussion and Conclusions 

The major task for present and future dark energy surveys is to determine whether dark en- 
ergy is evolving in time. Using the latest cosmological datasets (SN, CMB and LSS), we have 
performed a Bayesian analysis to extract the general form of the dark energy equation-of-state 
parameter, employing an optimal nodal reconstruction where w{z) is interpolated linearly 
between a set of nodes with varying Wz^-ysilues and redshifts. Our method has the advantage 
that the number and location of nodes are directly chosen by the Bayesian evidence. We have 
also explored standard parameterisations which include the CPL, JBP and FNT models. We 
find our results to be generally consistent with the cosmological constant scenario, however 
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Figure 5. Left: ID marginalised posterior distributions of the standard cosmological parameters, of 
each corresponding model listed in the right table. Right: comparison of the Bayes factor Bi^\ for 
some selected models with an extra- number of parameters Npar- Each description is compared respect 
to the ACDM model. 



the dark energy does seem to exhibit a temporal evolution, although very weak. Besides 
the cosmological constant, the preferred w{z) has w < —1 at the present time and a small 
bump located at z ~ 1.3, whereas at redshifts z > 1.5 the accuracy of current data is not 
enough to place effective constraints on different parameterisations. It is also interesting to 
note the presence of a narrow waist in many models, situated at z ~ 0.3, which is where 
the constraints on w{z) are tightest. A dominant feature throughout the reconstruction is 
the presence of the crossing of the PDL w = —1, obtained within the range < z < 0.5. 
Within the GR context, this phantom crossing cannot be produced by single (quintessence 
or phantom) scalar fields. Hence, if future surveys confirm its evidence, multiple fields or 
additional interactions should be taken into account to reproduce this important behaviour. 

All the models considered share a consistent set of primary cosmological parameters: 
f^b^^; ^DM^i^) 0, T, Hs, As, in addition to secondary parameters: Asz, Ap,Ac. The marginalised 
posterior distributions for these parameters are consistent with those obtained using only the 
concordance ACDM model. In Figure 5, we plot ID posterior distributions of the cosmological 
parameters for some selected models. We observe that their values remain well constrained 
despite the freedom in w{z). The only noticeable change is in the dark matter parameter, 
where the ACDM model displays the tightest constraints. In the same figure we include the 
corresponding Bayes factors, all of which are quoted relative to the cosmological constant 
model. The preferred Bayesian description of the w{z) is provided by the ACDM model, 
followed by the two-internal-node model Z2, introduced in this work. It is important to note 
that the CPL and JBP models, each with two parameters, are not able to provide an adequate 
description for the behaviour of w{z), and are hence strongly disfavoured using the priors 
chosen. The FNT model with four parameters, from which two of them remained uncon- 
strained, is significantly disfavoured. We stress that for the smallest prior range, the Bayes 
factor for the JBP model (which does not include the case wq = —1) is indistinguishable from 
that of the ACDM model, therefore pointing to a possible departure from the cosmological 
constant. 
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